Reconstructing the expansion history of the Universe with a one-fluid approach 
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Assuming that the Universe is filled by one single fluid, we present in the context of General Rela- 
tivity a possible explanation for the acceleration of the Universe. We use ordinary thermodynamics 
and the fact that small matter perturbations barely propagate in our Universe, to derive a general 
solution for a single fluid in which the speed of sound vanishes. We find a model that contains 
ACDM as a special case, and is compatible with current observational data. 
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To explain the (positive) acceleration of the Universe 
[H, it is usually assumed that, besides a dust-like fluid, 
the Universe is filled by an additional exotic fluid that 
accounts for about 75% of the energy density in the Uni- 
verse. Many models are known in the literature that in- 
tend to describe the nature of this additional component 
@; IH ! the simplest explanation is obtained by assuming 
the existence of a cosmological constant which is the ba- 
sic ingredient of the well-known ACDM model. Despite 
its simplicity, the ACDM model suffers from various the- 
oretical shortcomings Q, so it appears inadequate to be 
considered as a definitive model. In this work, we pro- 
pose a different approach based on the assumption that 
the dynamics of the Universe satisfies the laws of ordi- 
nary thermodynamics and on the observational fact that 
the speed of sound vanishes for a matter fluid Q. 

A perfectly homogeneous and isotropic cosmology with 
zero spatial curvature 1 is described by the Friedman- 
Robertson- Walker (FRW) line element 



ds 2 = dt 2 - a(t) 2 (dr 2 + 



r sin 



(1) 



Moreover, the gravitational source is assumed to be de- 
scribed by the energy-momentum tensor of a perfect fluid 
6]: T pv = diag(p(i), —p(t), —p(i), —p(t)). One can show 
that for a cosmological model satisfying the above sym- 
metry conditions it is possible to apply the laws of ordi- 
nary thermodynamics in a consistent manner Q • Let us 
consider a reversible and frictionless universe, satisfying 
the first law of thermodynamics, dQ — dll+pdV, and the 
equation of state (EoS) of an ideal gas p = p(C p — CV)T, 
where the heat capacities are defined as Cy — (fy) v 
and C p — , with h = U + pV being the enthalpy 

of the system. Furthermore, we assume that the evolu- 
tion of the Universe is adiabatic and reversible so that a 
polytropic relation holds, p — pp 1 Q, where 7 = C p /Cy. 
It is well known that currently the perturbations of the 



density, Sp, barely propagate in the Universe |5jj . It then 
follows that the speed of sound c 2 — can be con- 

sidered as vanishing. On the other hand, from the above 
thermodynamic assumptions we obtain for the Universe 
c s = y/j(C p — Cy)T so that for a vanishing sound speed 
the allowed solutions are formally C v — Cy = or 7 = 0. 
The first solution implies a pressureless universe whereas 
the second solution leads to a heat capacity C p = 0, 
equivalent to h — const, and a pressure p = —pCyT 
which is negative for positive values of the heat capacity 
Cv 

We conclude that the important result of assuming an 
ideal-gas-like universe with zero speed of sound is that 
the pressure in that universe can be zero or negative 2 . 
Moreover, let us notice that this is a consequence of the 
fact that each particle of a given fluid undergoes an early 
isentropic process, when the wave amplitude is infinitesi- 
mal. This condition follows from the fact that, in general, 
the entropy is proportional to the square of the velocity 
and temperature gradients. We will now consider a van- 
ishing speed of sound in the context of a FRW cosmol- 
ogy. Without putting any further information into the 
Einstein equations, we expect as a result a cosmological 
model in which the total pressure is negative and con- 
stant. Considering for the sake of simplicity the redshift 
z as the "time" variable defined by dz/dt = —(I + z)H, 
where t is the time coordinate, the conservation law for 
a generic fluid with EoS p = w(z)p, is 



dp 3(1 



w) 



dz 1 + z 
Assuming c 2 = dp/ dp = 0, we get 
dw 3w(l + w) 



p = 0. 



dz 



1 



z 



whose solutions are 

w = and w = — 



l-£(l + z) 3 



(2) 



(3) 



(4) 
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1 We limit ourselves to this case to be in accordance with observa- 



tions. The generalization to the case of nonzero spatial curvature 
is straightforward. 



2 Similar results were obtained by analyzing the ideal gas in the 
context of geometrothcrmodynamics 
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Notice that as a consequence of the thermodynamic hy- 
pothesis presented above, we obtain in a straightforward 
manner a dust-like term and an additional term with a 
time-dependent barotropic factor. This is a result of our 
model which is otherwise usually postulated arbitrarily 
in cosmology fl(ij |. 

Introducing the above solutions for w into Eq. @, we 
obtain the most general density for one fluid satisfying 
the condition c, = 



p{z) = { Pm + pi){\ + zf -p. 



(5) 



Notice that the term p m (l + z) 3 corresponds to the solu- 
tion w — 0. This general solution involves three constant 
parameters, namely, £, p m and p. 

Introducing p(z) into the first Friedmann equation 
H 2 = (8irG/3)p(z), we obtain the generic normalized 
Hubble rate, £=-#-, 

Ho 1 



X 



11 + 2 



fix- 



(6) 



where we adopt the convention fix = y- , with the critical 



density defined as p c = In Eq. 

condition E(z = 0) = 1, which gives 



l + fi 



A" 



fix 



|, we used the 



(7) 



where fl m = Therefore, to determine the Hubble 

Pc 

rate E(z) we need the constant fix, whereas to determine 
the barotropic factor w(z) in Eq. (j4j it is necessary to 
know the two constants fl m and fix- 

" 1 which follows 



Using the relationship as fix 



from Eq. (J7J), it is possible to rewrite Eq. as 



E= V«m(l 



fl A . 



(8) 



where SI, 
reads Oa 



and the cosmological constant term 



1 flm • 

The ACDM model is contained in Eq. ([5]) in the lim- 
iting case £ = with p = p m — 1 so that p m turns out 
to represent the sum of the baryonic and the cold dark 
matter densities. It follows that the model presented here 
is a generalization of ACDM and it is the result of the 
physical assumption that the Universe is made of only 
one matter fluid in which the the speed of sound is re- 
quired to vanish. Hence, the crucial difference lies in the 
fact that in the present model no cosmological constant 
is postulated a priori. 

In addition, difficulties related to the well-known prob- 
lems of coincidence and of fine tuning Q are solved in the 
context of the present model by the presence of a variable 
barotropic factor w(z) and by the fact that no cosmo- 
logical constant is assumed to be related to the vacuum 
energy. 

For instance, if we consider fl m = 0.274 and w — 
-0.980 [HI, we get from Eqs. Q and © £ w -0.02 



and fix ~ —0.712, which represents a value of fix, com- 
patible with an expanding Universe. 

Moreover, it is possible to infer the limits of the evo- 
lution of w(z) as 



w 



w 



|3J j for z = 0; 
= 0, for z 



(9) 



oo . 



Eqs. © show that at redshift z = u>o is a constant 
which predicts £ < for w > — 1. At higher redshift, the 
usual dust-like component dominates, because Woo = 0. 

If the Universe accelerates, the so-called acceleration 
parameter, defined as 



-1 - 



H 
IP 



(10) 



must be negative. From Eq. (jiop . for our model we 
obtain 



3 (1 + fi*) (1 + z) s 



2 + 2(l + f2 x )z 3 + z(3 + z) 



so that at z = it reduces to 



q = -(l + 3fi x ) 



(11) 

(12) 
-0.57 



For the particular value fix = —0.712, we get go 
which is in agreement with observations (12j |. 

At the moment in which the acceleration starts (q = 0), 
the correspondent redshift reads 



-1 



2{-fl x -2fi x -fl\) 



l + flx 

so that for fix = —0.712 we have z acc 



(13) 



0.7. 




FIG. 1: In this graphic is plotted q(z) for our model (dashed 
line) and ACDM (black line). The indicative values are Qm = 
0.274, fix = -0.712. 

Using Eqs. © and (|8|), it is possible to perform an 
experimental procedure to constrain the values of the 
constants fl m , fix and £. In particular, we employ the 
three most common fitting procedures: Supernovae la 
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(SNela), Baryonic Acoustic Oscillation (BAO) and Cos- 
mic Microwave Background (CMB). We will use of the 
most recent updated Union 2 compilation [ll|, which al- 
leviates the problem of systematics. 

Thus, associating to each Supernova modulus fj, the 
corresponding la error, denoted by cr M , we define the 
distance modulus fj, = 25 + 51og 10 j^:, where (Il(z) is 
the luminosity distance 



d L (z) = (1 + z) 



dz' 



/o H{z>)> 

and we minimize the chi square, defined as follows 

(/x theor _ ob.)2 



XSN 



E 



(14) 



(15) 





Q. m {BAO) 


n m (CMB) 


0.275 ±0.016 


0.264 ± 0.016 


0.280 ±0.019 


asN) 


t{BAO) 


£(CMB) 


0.036 ± 0.009 


-0.050 ±0.011 


-0.040 ±0.012 




Vt x {BAO) 


n x (CMB) 


-0.765 ± 0.058 


-0.700 ± 0.061 


-0.729 ± 0.050 



TABLE I: Summary of the numerical results for our model; 
the constants Q. m and £ have been fitted by using Eq. Q 
and the parameter Qx by using Eq. ([BJ. The chi square's 
are X SN.ia = 1-010, X Iao = 0.989, xcmb = 1-000, and 
XsNeia = 1-030, xbao = 1-010, Xcmb = 0-997, respectively 
for the test with Eq. (O and Eq. ©; the mean values for the 



three parameters are Q 



m (mean) 



0.273, Q 



X (mean) 



-0.731 



and ^mean = —0.042. Note that the value of Ho for the SNela 
tests is H = 72 ± 5 Kms' 1 Mpc~\ 



The second test that we perform is related to the ob- 
servations of large scale galaxy clusterings, which provide 
the signatures of the BAO [13| . We use the measurement 
of the peak of luminous red galaxies observed in Sloan 
Digital Sky Survey (SDSS), denoted by A 



A 



r h q - 


1 

3 


\ 1 / 


-H{z B Ao)- 




zbao Jo 



Hp 

H{z) 



dz 



(16) 

with zbao = 0.35. In addition, the observed A is esti- 
mated to be 



0.469 



0.95 
098 



-0.35 



(17) 



with an error a a = 0.017. In the case of the BAO mea- 
surement we minimize the chi square 



Xbao ~ 



A-A 

OA 



(18) 



Finally, for the CMB test we define the so-called CMB 
shift parameter 



R = 



H(z) 



dz . 



(19) 



with zcmb — 1091.36 pjj. It gives a complementary 
bound to the SNela data and BAO because the SNela 
redshift is z < 2, zbao — 0.35, while here z ~ 1100. 
We minimize the chi square 



Xcmb — 



R — R oi 

OR, 



(20) 



It is important to note that BAO and CMB do not 
depend on the values of Ho. We summarize the results 
of this numerical procedure in Tab. I 

The values of fi m , fix and £ of Tab. I are in agreement 
with the theoretical results showed previously. 

During the last decades, different parametrizations of 
w(z) were proposed @- For instance, w — w\ ± W2Z : 



w = w a + wp log(l ± z) or w = Wo + w a (l — a). In 
particular, the third case was introduced by Chevallier, 
Polarski and Linder and it is referred to as the CPL 
parametrization [l5l |. The CPL parametrization has the 
advantages that at low and very high redshift it reduces 
to constant values, respectively w(z — > 0) = wq, and 
w(z — > oo ) = wo + w a - However, all the parametrizations 
suggested so far for w(z) are either ad hoc proposals or 
the result of phcnomcnological assumptions only. Our 
model predicts a theoretical barotropic factor w(z), and 
it is also able to reproduce previous results (see, for in- 
stance, (U). 

Moreover, the barotropic factor is also connected to an 
interesting quantity, that one might consider as a natu- 



ral measure of time variation, namely, 



dw 



dln(l+z) 



111 



models involving a scalar held tp with potential V((p), this 

quantity is related to the slow- roll potential cx -y- , in the 
region 2 = 1, where the scalar field is most likely to be 
evolving as the epoch of matter domination changes over 
to dark energy 3 domination. For the CPL parametriza- 

= ^f-, while for the present model it 



tion 



dw 



d\n(l+z) 



dw 



d\n(l ± z) 



-24ft 



l + O 



A' 



z=l 



X " 



7n 



X 



(21) 



By comparing our result with CPL we have two solutions. 
One of these solutions is physically compatible with the 
accelerating scenario; in fact, by considering £l m (mean) 



15| . we find fix 



and the indicative value w a = 0.58 
—0.716, in agreement with the observational results. 

We present below the graphics of the evolution of w(z), 



3 It is a common opinion to refer to the missing ingredient, driving 
the acceleration, as dark energy. 
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q(z) and the expansion history for a(t), given by (l5l Il6j 



-0.8 




-1.0 ' 

2 4 6 8 10 12 14 

FIG. 2: In this graphic is plotted w(z) (Y axis) for our model 
(dashed line) and CPL (black line). The indicative values are 
n m = 0.274, fix = -0.712, too = -0.93 and w a = 0.58. 
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The presence of baryonic and dark matter is generally 
intertwined with the addition of an exotic fluid which 
drives the acceleration. Unfortunately, all the attempts 
to describe this unexpected acceleration, suffer from var- 
ious shortcomings. Moreover, the ACDM remains the 
favorite fitting model to describe the Universe dynamics, 
by including in Einstein equations a second fluid char- 
acterized by the cosmological constant. We showed that 
it is possible to discard the existence of a second fluid 
and, instead, to use ordinary thermodynamics. Assum- 
ing a vanishing speed of sound in order to guarantee that 
small matter perturbations do not propagate, we obtain 
a theoretical parametrization for w, which generalizes the 
ACDM model, reducing to it in a special case. In addi- 
tion, the model presented here solves the coincidence and 
fine tuning problems in a straightforward manner. Using 
the cosmological tests of SNela, BAO and CMB, it was 
shown that our model is able to reproduce the observable 
Universe, and is in agreement with the theoretical limits. 
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